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ABSTRACT 


UPDATED POSITION AND ICE VELOCITY 
FOR THE AIDJEX MANNED CAMPS 


A.S. Thorndike and T.O. Manley 


During 1975 and 1976, measurements of sea ice motion were made as 
part of the Arctic Ice Dynamics Joint Experiment (AIDJEX). In an earlier 
report (Thorndike and Cheung, 1977), energy at the inertial period was at- 
tenuated by approximately 50% due to the choice of a low covariance factor, 
Q = 100, during Kalman filtering. 

For oceanographic investigations, the data set proved to be unac- 
ceptable particularly in the summer months when there is a significant amount 
of energy at the inertial period. 

This report is the result of the reprocessing of the original position 
data taken from the four manned camps (Caribou, Blue Fox, Snowbird and 
Big Bear) with a covariance factor Q = 100000. 

The higher Q allows data of the inertial period and longer to pass 
through the filter with less than 1 percent attenuation. 

Because of the large amount of information contained in the data set, 
it was divided into two parts. This report (Volume 2) provides position and 
ice velocity from 18 October 1975 to 4 May 1976. Volume 1 provides position 


and ice velocity from 11 April 1975 to 17 October 1975 (Thorndike and Manley, 


1980). 
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Introduction 


The objective of the Arctic Ice Dynamics Joint Experiment, as set 
forth in its most general statement, was to reach an understanding of the 
dynamic and thermodynamic interaction between the ice cover and its en- 
vironment (Maykut et al., 1972). One component of the experiment was an 
array of drifting stations at which frequent measurements of the position 
were made (Figure 1). For the present report those measurements have 
been edited using Kalman smoothing techniques to give the best estimates of 
position and velocity. (The raw measurements themselves are not reported 
here, but are available through the Polar Science Center.) The data set 
described here is restricted to the main AIDJEX experiment, which began in 
March 1975 and ended in May 1976. 

This data report represents the most complete documentation of ice 
motion available, but it is limited to one year and to the Beaufort Sea. Par- 
ties interested in earlier data, more widely dispersed in space and time, 
should contact the Polar Science Center, University of Washington, Seattle, 
Washington, which has compiled most of the ice drift data beginning with the 


drift of the FRAM in 1893 to the present. 
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Data Acquisition 


A full description of the measurement systems used in this program 
is beyond the scope of this report. The list of references includes sources 
of more detailed information (Wasilew and Vivian, 1976; Martin, P., 1978). 
All of the hardware we used was supplied by the Satellite Positioning Cor- 
poration, Houston, Texas. 

Four manned camps -- identified by radio call names Big Bear, Cari- 
bou, Blue Fox and Snowbird -- were tracked using NavSat (Navy Navigation 
Satellite System). ‘The geographical position of each of the four camps was 
measured at unevenly spaced times using Doppler satellite navigation tech- 
niques. The raw measurement is the frequency of a signal transmitted from 
the satellite to the ice station. The measured frequency is affected by a 
Doppler shift related to the rate of change in distance between the satellite 
and the ice station, which is itself related to the unknown ice coordinates 
and the known satellite orbit. During each satellite orbit several frequency 
measurements were made and from them an overdetermined solution was found 
for the unknown coordinates. Data recording in the field was on magnetic 
tape with hard copy printout for backup. 

A constellation of six polar orbiting satellites was maintained for this 
system. The orbital period is about 107 minutes. An observer at high lati- 
tudes can track every orbit. Of the 81 potential passes each day, many are 
unusable (because of radio interference between overlapping passes, for in- 
stance). We averaged about 25-30 usable passes per day at each station. 


Single channel (400 MHz) receivers were used to count the Doppler cycles 


for roughly 30-second intervals. The duration of each Doppler counting 
interval was measured independently. The frequency standard used in the 
receivers was stable during a 20-minute pass to about 1 part in i. 

Fix solutions were done in real time at each manned camp. The 
solution used only the Doppler data from the 5 minutes on either side of the 
time of closest approach of the satellite. No additional editing was performed. 
A correction to each Doppler count was based on the duration of the counting 
interval to account for the random variations in the counting interval. The 
correction also included the fixed timing delay in each receiver. The fix 
solution ignored the effect of the velocity of the ice camps. 

In fact, two receivers with antennas separated by about 100 meters 
were operated at each manned camp. In addition to serving as a backup, 
the measurements at the second antenna were used to define the azimuth of 
the line joining the two antennas. This was possible when both receivers 
produced acceptable fixes from the same satellite pass. Estimates of the 


azimuths obtained at each of the manned camps were reported by Thorndike 


and Cheung, 1977. 


Preprocessing 


The original data tapes from the manned stations measurement sys- 
tems typically contained several "glitches" per tape because of power failures 
at the camp, or operator error, or some other reason. These tapes were 
transcribed discarding those records which had either a parity error or an 
incorrect number of words. The tapes which resulted from this process 
were read without error on the CDC 6400 at the University of Washington, 
where the bulk of the processing was done. 

Three problems had to be dealt with to prepare the sequence of 
measurements at each station for the actual Kalman filtering. First, data 
recording problems had created occasional gaps in the manned camp data. 
At times the backup teleprinter records from the field could be used to fill 
these gaps. This was done by keypunching the fixes directly from the 
printed records. Second, the data sequence had occasional points which 
appeared to be out of sequence in time. This was sometimes correctable, 
as when the operator had keyed in the wrong date, but usually the situation 
was ambiguous and points out of chronological order were dropped. The 
third and most serious problem involved bad fixes. A glance at a plot of 
the raw time series (Figure 2) illustrates the problem. Because outlying 
fixes have a deleterious effect on the Kalman filtering results, an attempt 
was made to remove them during preprocessing. The algorithm employed 
compares each fix in latitude with the median of the latitude of the 10 fixes 
preceding and the 10 fixes following it. When the differences exceeded a 


preset tolerance of 10 kilometers, the fix was eliminated. The algorithm 
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was applied twice, with a smaller tolerance of 1 kilometer the second time. 
The same procedure was then applied for longitude. 

Visual checks were an important part of the procedure. Plots of 
the data points were produced before and after preprocessing. At times it 
seemed from the plots that we were not detecting certain bad fixes, or that 
we were throwing out good ones, and it was necessary to adjust the toler- 
ance levels. Edited raw data plots like the one shown in Figure 3 provided 
assurance that the data were ready for filtering. 

For the main experiment, AIDJEX adopted the convention of number- 
ing days consecutively, beginning with day 1 = 1 January 1975 and ending 
with day 500 = 14 May 1976. Throughout this data report, this numbering 
system is used quite extensively. Appendix 3 provides a conversion table 


from AIDJEX days to calendar days. 
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Kalman Smoothing 


In the smoothing used here an assumption is made regarding the 
motion of the ice. If the state of the ice at time in is represented by a 
vector Xn , containing, say, the position, velocity, and acceleration of a 
point, then the state at some future time is assumed to be partially deter- 


mined by Xn : 


X = POSS OHA | Bari (1) 


where © and I are matricies known from the physics of the situation, and 
be +1 represents an unknown random perturbation. In Kalman's formulation, 


which we follow, the random effects are assumed to be Gaussian and White: 


Q ifn=m, and 
cov (W_W_) = (2) 
sein 0ifn#mMm 


We have chosen @ = 100000 m“hr °. 


Qualitatively choosing a small value of 
Q is equivalent to assuming that the ice moves in a completely predetermined 


way. (See from equations 1 and 2 that Q = 0 would imply that x, = 6 Xn 


+1 
a large value of Q corresponds to assuming that the ice experiences large 
sudden changes in its acceleration.) Thorndike and Cheung (1977) used 

Q = 100 m?hr ©, The lower value of Q assumed smooth behavior, but was 
unnecessarily conservative. The present value of Q was determined by trial 
and error. A sequence of raw fixes was used as input to the Kalman filter. 
Several different runs were done using different values of Q. The present 


value of Q appeared to give the best agreement between the output from the 


filter and the raw fixes. 


The original Q (100 mhr ®) used by Thorndike and Cheung removed 


approximately 50% of the inertial period amplitude as can be seen in Figure 4. 
A constant amplitude sinusoidal function whose period varied linearly with 
time was chosen as the input to the Kalman filter with Q = 100 m*nr ©, In- 
put periods ranged from 0.5 hours (at day 10) to 60.4 hours (at day 15) with 
the inertial period being located at day 10.9. An approximate 50% attenuation 
of the input function at the inertial period can be roughly estimated from the 
graph. 

Figure 5 shows the same input function and the resulting Kalman 
output with Q chosen to be 100000 m7 hr ®, The inertial period passed through 
the filter with very little attenuation. Another test was then run, this time a 
constant amplitude, constant period (12 hours) sine wave was used to simulate 
the inertial period of the Arctic latitudes. The results indicated that there 


was only a 0.03% attenuation from the original input function. 


For the Kalman reduction, we defined A = ae = ty and the 


i 
matrices 
x = position O: gg 7h Lah a2/2 
velocity ea ye (3) 
acceleration 0 O 0 
r = |a%e6 
a7 /2 (4) 
A 
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An identical equation can be written for the y-coordinates. The 
two coordinates are treated independently. 


The measurement Zn is related to X by 


Zn = (1 0 OFX Er 


: s : 2 
where VES is a measurement error with assumed covariance R = (60 meters) . 


The response of the measurement and processing system was investi- 
gated as follows. Four hundred sample points were drawn independently 
from a Gaussian distribution with zero mean and variance (50 meters) °. A 
sampling interval of one hour was assumed. This noise simulated the Nav- 
Sat measurements. The sequence was used as input to the Kalman filter 
procedure. The velocity output from the filter shows two effects. First, 
since velocity is the time derivative of position, we expect to find the spec- 
trum Sw) behaving as Ms S,.(w) , where 5, () is the spectrum of position. 
S, (w) was assumed to be white; S,. (w) = a , which implies S,(w) = a o : 
The second effect is the tendency of the filter to smooth out high frequen- 
cies. Qualitatively then, we have something like the diagram in Figure 6. 

The actual spectrum of ice velocity is shown in Figure 7. The dia- 
gram illustrates how the measurement errors and finite sampling rate limit 
the accuracy of the velocity estimates. Most of the variance in the velocity 
errors occurs for the frequencies greater than 3 cycles per day. For lower 


frequencies we can estimate the contributions to the variance as 0 - 1.5 cycles 
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Figure 6 - Diagram showing qualitatively, the response of 
the measurement and processing system. 
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Figure 7 - Spectrum of ice velocity from a test run on the 
processing system. For a more complete de- 
scription, see text. 
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per day: o* = 3x10 


2 


6 2 


m? sec“, 1.5- 3.0 cycles per day: o© = 2x10 


total variance 0 - 12 cycles per day: o” = iit) = m2 pechie 


5 
m” sec” 
The general smoothing problem is to estimate X at time th using the 
entire set of observations (Z;) i= 1....m, and to give the error variance 
of that estimate. The solution equations, due to Kalman, are reproduced 
in Appendix 1. All positions were converted to an X-Y coordinate system 
prior to actual filtering. This coordinate system along with its governing 
equations are shown in Appendix 2. Typical estimation error variances for 


position and velocity are 40 meters and 1.0 cm/sec respectively. 
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Results 


The final data listing consists of five major columns. The first 
column is the time in decimal AIDJEX days corresponding to the observation 
of the position and ice velocity at a particular camp. The remaining columns 
indicate the camps for which data are applicable. Subdivisions under each 
major camp heading are LAT (latitude in decimal degrees), LONG (longitude 
in decimal degrees), and the north and east components of ice velocity in 
Cm/sec, NIV and EIV respectively. 

Due to the large amount of information contained in the data set, it 
was divided into two parts. This report (Volume 2) covers AIDJEX days 
291 (18 October 1975) to day 490 (4 May 1976). Volume 1 covers AIDJEX 
days 101 (11 April 1975) to 290 (17 October 1975) (Thorndike and Manley, 
1980). | 

Figures 8- 11 show a more detailed drift track of each of the manned 
camps. An insert shows where the region lies with respect to the north slope 
of Alaska and Banks Island. Asterisks superimposed on the drift track indi- 
cate the positions at integral multiples of 20 days. The beginning and end- 
ing days are noted for each trajectory. A dashed line indicates a long period 
of missing data. The plotted region is 500 x 500 km and is aligned with the x, 
y coordinate system shown in Appendix 2. Figures 8 - 11 were reproduced 
from Thorndike and Cheung (1977). 

Table 1 lists the last known positions and ice velocities from camp 
Big Bear during the spring of 1976. Big Bear was the original main camp 


but was abandoned in early October of 1975 due to ice breakup. In spring 


16 


of 1976, RAMS buoy number 604 was placed on the site of the camp which by 
this time had been spread over a few kilometers. Data cover from day 465 
to 485. Because of RAMS buoy lower sampling rate and higher noise prob- 


lem, the data in Table 1 were not reprocessed. 
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APPENDIX 1 


KALMAN SMOOTHING EQUATIONS 


Symbols 


x = state vector (n x 1) 
d = transition matrix (n xn) 


T = coefficient matrix for random perturbations to physical 
model (nx m) 


= white noise driving the physical process (m x 1) 

= covariance matrix for W (m x m) 

= measurement vector (2 x 1) 

matrix relating measurement to state vector (2% x n) 
= measurement error vector (2x 1) 


= measurement error covariance (2 x 2) 


vad otnod = 
iI 


= prediction, filtering and smoothing error covariance 
matrix (n x n) 


= identity matrix (n xn) 

= Kalman gain matrix for filtering (n x 2) 
smoothing gain matrix (n x n) 

= number of measurements in the data sequence 
= index for the data sequence 


= time 


Oe ea ae 
i 


The superscript T denotes the transpose operation. The symbol 
XG lj) means the estimate of x at time t given observations (Z,, Zo» Rai z,) A 
The notation here follows Meditch (1969). 


Prediction 


Given observations (Z,, wort; Za); predict the state at the future 


time th +1? and find the prediction error covariance. 


rr a a a ee Ped 
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X (k + 1, | k) o(k+1, k) X(k|k), A(1) 


P(k+1|/k) = (k+41,k)P(k|k) 0? (k+1,k) 
+ T(k+1,k)Q(k) IP? (k = 1,k). A(2) 
Filtering 
Given observations (Z,; eer 2.) and a new observation Zac+1)? 
estimate the present state at time ti. +1 and find the filtering error co- 
variance. 
XUk+1]/k+1) = X(k+1]k)+K (K+ 1) (ZK+D 
- H(k + 1) x(k+1|k)), A(3) 
P(k+1|[k+1) = (1- K(k + 1)H(k+1)) P(kK+1]{k), A(4) 
where 


K(k +1) = P(k+1/k) H) (k+ 1) (H(k + 1) P(k + 1]k) H! (k + 1) 
1 


+ R(k+1))— A(5) 
Smoothing 

Given the complete data set (Z,; eens Zy)> estimate the state at 
any point the and find the smoothing error covariance. 
K(k |N) = Xk k) +A Ck) (X(K + 1[N) - X (K+ 1|k)) ‘A (6) 
P(K|N) = P(k|k + A(k) (P(k + 1|N) 

- P(k +1]k)) A? (k), where A(7) 
Aka P(k|k) 7 (k + 1,k)P b(k+1|k). A(8) 
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Implementation of the Filter Algorithm 

Starting from an initial condition x (0|0), P(0|0), a prediction is 
made to time t,. The prediction error covariance matrix P (1| 0) is used to 
obtain the Kalman gain K(1). The filtered estimate X(1| 1) is then found 
by adding to the predicted x (1| 0) a quantity which is proportional to the 
difference between the observation Z(1) and H(1)x(1|0). H(1)x(1|0) is 
the value the measurement would have if the true state were actually X(1| 0) 
and no measurement errors were made. The algorithm is repeated as each 
new data point is entered until the last filtered quantities x (N | N) and 
P(N |N) are produced. 

The smoothing algorithm, proceeding backward in time, utilizes 
results from the prediction and filtering calculations. The first step is to 
find X (N-1]N) and P(N-1|N) from equations A(6), A(7) and A(8). No 
direct use of the measurements is made in smoothing. The algorithm simply 
corrects the filtered vector x (k|k) by a quantity proportional to the differ- 
ence between the prediction vector at time ted and the smoothed value at 
time tk +1° 

In the calculations done for this application, the CDC 6400 computer 
time required to do the prediction, filtering, and smoothing is about 0.5 sec- 
ond per data point, nearly half of which is absorbed in the matrix inversion 
in A(8). Because of the recursive nature of the filter, demands on central 


memory storage space are very modest. 


Reproduced from Thorndike, A.S., 1973, An integrated system for meas- 
uring sea ice motions. In Ocean 73, IEEE International Conference on En- 
gineering in the Ocean Environments, IEEE publ. no. 73 CHO 774-0 OEC, 
pp. 490-499. 
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APPENDIX 2 


COORDINATE SYSTEMS 


(km) 
- 1500 -1000 


-180° 


-165° 


Position measurements were made in geographical coordinates 
(latitude north, longitude east). The smoothing operation was 


done in a Cartesian system (x,y), where 


* 
I 


110.949 (90° - latitude) cos (longitude) (km) 
110.949 (90° - latitude) sin (longitude) (km) 


td 
TT 
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APPENDIX 3 


CONVERSION TABLE FOR AIDJEX DAYS 
TO CALENDAR DAYS 


For the main experiment, AIDJEX adopted a convention of number- 
ing days consecutively, beginning with day 1 = 01 January, 1975 and ending 
with day 500 = 14 May, 1976. 

In the conversion table, the first column is the AIDJEX day, the 
second is the corresponding day of 1975 or 1976 and the third entry is the 


calendar date. 
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Appendix 3 - Conversion Table 
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TABLE 1 


Last known positions and ice velocities for camp Big Bear. Taken 
from RAMS buoy R604 (Thorndike and Cheung, 1977). Days are in decimal 
AIDJEX days. Negative longitudes imply west longitude. NIV and EIV are 
north and east ice velocities respectively. Negative values in the NIV or 


EIV columns imply south or west velocities. 
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TABLE 1 


Last known positions and ice velocities for camp Big Bear. 
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TABLE 1 (cont.) 


Latitude Longitude 
72.9226 -141.3516 
72.9166 -141.2966 
72.9143 —141 2500 
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72.9459 =La Li T2 
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12.9188 -141. 3874 
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TABLE 1 (cont.) 


Latitude Longitude 
72.9276 -141.4327 
72.9388 -141. 4558 
72.9441 -141.4542 
72.9443 -141.4557 
72.9406 -141.4666 
72.9345 -141.4747 
72.9255 -141.4718 
72.9170 -141.4579 
172.9123 -141.4397 
12.9125 -141.4319 
72.9174 -141.4436 
72.9195 -141.4505 
72.9192 -141.4565 
72.9183 -141.4761 
72.9191 -141.5241 
72.9218 -141.5924 
72.9266 -141.6676 
72.9324 -141.7305 
72 ..9393 -141. 7560 
72.9489 =1417,7819 
72.9674 -141.8250 
72.9790 -141.8889 
72.9965 -141. 9578 
73.0138 -142.0193 
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TABLE 1 (cont.) 


Latitude Lon gitude 
fo.UoLU -142.0579 
73.0449 -142.0749 
73.0496 -142.0779 
73.0499 -142.0832 
73.0506 -142.0979 
73.0543 =142.1235 
73.0561 =142.1321 
73.0542 -142.1247 
73.0487 -142.0968 
73.0373 -142.0555 
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Position and Ice Velocity Data 


The following section comprises the updated positions and ice veloci- 
ties for the four manned camps Caribou, Blue Fox, Snowbird and Big Bear. 


Data span from 18 October 1975 to 4 May 1976. 
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